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Abstract 

If neutrino mass is a function of the Higgs potential then minimum of the total thermodynamic 
potential (which is the Higgs potential minus the neutrino pressure) can shift from the standard 
electro-weak vev v = 246.2 GeV by a small amount which depends on the neutrino pressure. If 
the neutrino mass is a very steep function of the Higgs field then the equilibrium thermodynamic 
potential can act like the dark energy with ~ — 1. Choosing the neutrino mass as logarithmic 
function of the Higgs field and a heavy mass scale, we find that the correct magnitude of the 
cosmological density of the present universe px ~ (0.002ey)^ is obtained by choosing the heavy 
mass at the GUT scale. 
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1. INTRODUCTION 



The numerical coincidence in the cosmological dark energy density px = (0.002eV)^ and 
the neutrino mass scale Am^ = 2.5 x lO^^eV^ which solves the atmospheric neutrino puzzle 
has led to the ideas that the minimum of the dark energy potential is set by the dark 
energy dependent neutrino mass 1-3]. There is a numerical relation between the GUT 
scale M = lO^^GeV the weak scale v = 246.1GeV and the dark energy density px ~ v^/M"^. 
This motivates us to use relate the Higgs potential to the dark energy. Since neutrinos are 
the most abundant particles with mass, it is natural to use the neutrino-Higgs interaction 
to shift the Higgs field from its minimum of the potential = (A/4)(0^ — v"^). The pressure 
of the neutrino fluid depends on the neutrino mass. If the neutrino mass is a function of the 
Higgs field then the field rolls down to the lowest point in the total thermodynamic potential 
i.e the Higgs field dynamically adjusts to a value which maximizes the total pressure of the 
Higgs-neutrino fluid. Maximizing the total pressure results in the thermodynamic potential 
minima shifted from (0) = v to (0) = v + am- For this scenario to work and give the correct 
dark energy density the neutrino mass as function of the Higgs field must have a very steep 
gradient. We have considered logarithmic neutrino mass functions which meet this criterion. 
The neutrino mass function has a mass scale M, and dimensionless small parameter e as free 
parameters. To give the correct dark energy density it turns out that M must be close to the 
GUT scale owing to the numerical relation (pxY^^ ~ v'^/M. Neutrino mass in this scenario 
turns out to be = 2 x 10~^eV if neutrinos are assumed to be non-relativist ic. This is the 
range consistent with the solution of the atmospheric neutrino problem. One well known 
problem of the neutrino couplings with low mass quintessence type fields is that they give 
rise to a long range attractive force between neutrinos which could lead to an instability as 
pointed out in reference I?]. In our scenario the Higgs coupling will give a contact interaction 
between neutrinos which does not lead to formation of neutrino nuggets and clustering at 
astrophysical scales. 



2. THERMODYNAMIC POTENTIAL OF HIGGS-NEUTRINO SYSTEM 

The Higgs potential is assumed to be at zero temperature and the thermodynamic poten- 
tial equilibrium between the Higgs and the neutrino pressures is due to the Higgs dependence 



2 



of the neutrino mass which affects the neutrino pressure. In the absence of the neutrino cou- 
phng the Higgs field attains a vev ($) = vj \/2 such that the Higgs pressure = 0, 

The Higgs field $ is assumed to be at zero temperature and its thermodynamic potential, 
in the unitary gauge where $ = ^(0, f + 0)^, is 

n, = -p, = ^i<i>'-vr (1) 

where v = 246.2 Gev is the scale of electroweak symmetry breaking. When the neutrino 
mass is a function of the Higgs field then the net thermodynamic potential 

n = n,{(t>) + (2) 

can have a minima at a value of (pm 7^ v. If m({/)) is of such a form that at the minima the 
thermodynamic potential is positive, > then we have an effective negative pressure 

from the Higgs-neutrino interaction. The thermodynamic potential has a minima at = 0^ 
where 

= n^{(f)j' + m{(f)^y = (3) 

where prime denotes derivative w.r.t 0. We look for a functions m(0) for which the solution 
(pm of equation ([3]) is such that 

^<t>{(J)m) + ^u{(l)m) > (4) 

We expand the Higgs field close to the minima of Q^, by taking (f) = v + a. If the minima 
occurs at a non-zero am then we know that Qfj,{am) > (from the form of the Higgs 
potential). We also know that the pressure of the neutrino fluid is always positive which 
means that Q^{am) = —Pu < 0. Moreover the magnitude of the pressure Py of relativistic 
neutrinos is five orders of magnitude smaller than the dark energy density p\ = (0.002eV)^ 
(and even smaller for non-relativistic neutrinos). This implies that at the minima of Q the 
Higgs potential must dominate over the neutrino potential fl^{am) = ^'v'^cTm ^ ^u{o'm) 
therefore 

Pa = At;V^ (5) 
implying that the Higgs minima should be shifted from = f by the amount 

CT^ = ±^1.6 X 10"^ W (6) 
V A 
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in order to explain the observed cosmological acceleration. The neutrino mass function 
m(0) must be a steep function of the the Higgs field, m(crm)' ^ 1 in order to satisfy the 
simultaneous requirements Q^{am,) ^ — ^i/(o"m) and fi(^((Tm)' + Qt^^am)' = 0. 
The effective equation of state of the Higgs-neutrino fluid is 

) 

as Pu,Pu ^ V(f>{o'm)- The Higgs mass = V^{am) — 4Af^ is much larger than the Hubble 
rate and the Higgs field stays at the bottom of the potential as the universe expands and 
the value of the minima a^n changes with cosmic time. The change in the pressure of the 
Higgs-neutrino fluid with time is related to the change in the am with time as 

_=(^_(_V,M + P„M)j__ (8) 

The first bracket in the rhs of equation ([8]) dynamically adjusts to zero at thermodynamic 
equilibrium. The parameters in the potential must allow the solution of (— + P^)' = in 
the present era. This explains the 'why now' question about why the dark energy dominates 
in the present era. Once the pressure dynamically adjusts to be stationary the density also 
is stationary owing as the equation of state w ~ — 1, 

dp dP „ . , 

-T - — r = (9) 

dt dt ^ ^ 

We note that the thermodynamic principle we apply is that the system adjusts to maximize 
the pressure (or minimize the thermodynamic potential ). A justification of this from first 



principles can be found in references 16Nl8|. We assume that the Higgs and neutrino 



pressures are coupled through the neutrino mass and the Higgs field in the potential adjusts 
such as to be at the minima of the thermodynamic potential Q = — P^- This principle 
of minimizing the thermodynamic potential is also advocated for the case of mass varying 
neutrinos in reference [lOj]. The thermodynamics of stable structures formed from mass 
varying particles is discussed in j2oj |. 

Neutrinos decouple from thermal equilibrium with the radiation fluid at a temperature 
T = IMeV when their distribution function is the ultra-relativistic Fermi-Dirac distribution, 
f{p) = (6xp(p/T) + 1)^^. Below the temperature of 1 MeV both temperature and momenta 
red-shift with the scale factor as 1/a and the distribution retains the same form even though 
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they at the present temperature they may be non-relativistic. The thermodynamic potential 
of neutrinos which decoupled when they were ultra-relativistic can be written as 

n„W = -P. = -j./|^i(^)/(p) (10) 

where E = \Jm[(j))'^ + p^. Here /(p) is the statistical distribution function of neutri- 
nos. If neutrinos have only weak interactions with other particles and decouple from 
thermal equilibrium at T ~ IMeV then their distribution function remains relativistic 
f{p) = (6xp(p/T) + 1)^^ even in the present era. One can evaluate the thermodynamic 
pressure of neutrinos whether they are relativistic or non-relativistic in the present epoch 
by evaluating (fTOl) in the relativistic E = p + m? / {2p) or non-relativistic E = m + p^ / (2m) 
limits. 

Beacom et al 



21| have raised the possibility that if neutrinos have a large interaction 
with some hght boson then they do not decouple at T ~ IMeV but can continue to be 
at thermal equilibrium with the radiation bath of the bosons till the present epoch. This 
possibility is not ruled out by big-bang-nucleosynthesis and CMB observations. The neutrino 
distribution function in this case would be f{p) = {exp{E/T) + 1)"^ and the pressure of the 
thermally coupled neutrinos can be evaluated in the for the two possibilities that neutrinos 
are either relativistic or non-relativistic in the present epoch. We study the case of un- 
decoupled neutrinos as an interesting example of situation where the neutrino pressure is 
different compared to the decoupled neutrinos (although the decoupled neutrinos also have 
a relativistic-thermal distribution) . In the following sections we give examples of m(0) for 
the case of relativistic and non-relativistic neutrinos which can give rise to a net negative 
pressure of the required magnitude. 

3. NON-RELATIVISTIC NEUTRINOS 

We choose the neutrino mass function as 

m(a) = -^\og[^{l-.)') (11) 
where M is some large mass scale, v = 246.2GeV and e is a small dimensionless parameter. 
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3.1. Non-relativistic decoupled neutrinos 



= -Pu = -g^ j 



r £p 1 


/p2 p4 X 
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/ (27r)3 3 ' 
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' exp(p/r) + 1 



45C(5)T5 2835 C(7) 

+ o \ (12) 



27r^ m Svr^ 

where we have taken the distribution function of neutrinos to be relativistic, and we have 
taken Qi = Q assuming that all three species of neutrinos are non-relativistic in the present 
epoch. 

Then the total thermodynamic potential (to the leading order in m/T) is 

nw = A„V-iiiMl! (13) 

^ ^ 27r2 m ^ ^ 

has minima at the value of cr = am given by the equation Q' = 0, 

2Xv^am = ^ : — 772 r (14) 

This equation can be re-written in terms of the neutrino mass m{am) (equation ( |TT|) ) and 
the dark energy density px (equation ([5]) ) as, 



( - .)) . -liiiai! 1 . _,o X 10- eV (15) 

A solution of this equation can be obtained (in the approximation e ^ 1), by choosing 

= 1' (16) 

and using log(l + x) ~ x for small x, to get 

m(a^)e = 1.0 X lO^^eV (17) 
Writing ( !T6ll in terms of the dark energy density ([5]) we have, 

PA = ^ (18) 

from f|T8l) we see that taking M = 3 x lO^^GeVA^^^ we get the correct magnitude p\ = 
(0.002eV)^ for the dark energy density. The neutrino mass (fTTl) is 

, , /3 X IQi^GeVAi/n , , 

m(a„) = - = 0.002eV I J (19) 



The two observable quantities p\ and m depend on only one parameter M. The other 
parameter determines the temperature at which the VL minimization takes place. If m = 
0.002eV then from the minimization condition fll4p and (11 7p we see that taking e = 0.5 x 10^^ 
ensures that the minima of the thermodynamic potential occurs at the present era (when 
the neutrino temperature is T = 1.7 x 10~^eV). 



3.2. Non-relativistic non-decoupled neutrinos 

If neutrinos have a large interaction with some light boson then they do not decouple at 
T ~ IMev but can continue to be at thermal equilibrium with the radiation_^bath of the 
bosons till the present epoch. This possibility was studied by Beacom et al [21] who found 
that this possibility is not ruled out by big-bang-nucleosynthesis and CMB observations. As- 
suming that the bosons- neutrino fluid was relativistic at some temperature T = (A/liy^^T^, 
and subsequently the neutrinos are non-relativistic then from entropy conservation it follows 
that that the present common temperature of the non-relativistic neutrinos and relativistic 
bosons is T = = 3.8 x 10"^ eV Q- 

If the neutrinos are non-relativistic at the present epoch then the neutrino pressure at 
thermal equilibrium is describe by 

-n„ = p„ = 6(!!M)'%='^exp(z!!M) (20) 

The minima of the thermodynamic potential occurs at (fi^ + Q^)' = 0, and is given by the 
solution of 



- ^m(a„)V^T^/^ ( ^ _ !!^^ exp ( ( ^) = 
(27r)3/2 ^ T J ^ \ T J \eMamJ 

We choose the large mass scale M in the neutrino mass function ( |TT|) such that 
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(21) 



(22) 



which implies that M = 3.0 x lO^^GeVX^^^. Assuming the parameter e <^ 1 the expression 
( ITTj) then gives us the neutrino mass, 

m(a^) = ^ = ^ 2 X lO-^eV (23) 

With this value of neutrino mass the ratio m(am)/T = 5.2\^^/^ which is consistent with 
our assumption that neutrinos are non-relativistic. One can solve thermodynamic potential 



minimization equation ( 12T]) for the remaining parameter e to obtain e = 5.9 x 10 ^, which 
ensures that the thermodynamic minima is achieved in the present epoch (when the neutrino 
temperature is = (25/11)^/3 2^^ ^ g g ^ iQ-^eV). 



4. RELATIVISTIC NEUTRINOS 

Assume that at-least one of the neutrino masses is smaller than the neutrino temperature 
T = (4/11)^/3To = 1.7 X 10~^ eV. We consider two possibilities (a) neutrinos decouple at 
T IMeV and have a relativistic thermal distribution (b) they are in thermal equilibrium 



with the radiation in the present epoch as studied in reference [2l|]. For both the relativistic 
neutrino species we consider the mass dependence on the standard model Higgs of the form, 

where again M is some large mass scale, v = 246.2GeV and e is a small dimensionless 
parameter. 

4.1. Relativistic decoupled neutrinos 

The thermodynamic potential of relativistic neutrinos in terms of neutrino mass m((/)) is 
obtain from (fTOj) in the relativistic limit m < T is 

= -P. --'J'j§f\{"-%) e.p(p/V) + 1 

+ (25) 



360 72 

where we have taken Qi = 2 for one relativistic neutrino species. 
Then the total thermodynamic potential is 

^ ' 72 ^ ^ 360 ^ ^ 

The minima of fl occurs at the value of a = given by the equation Q' = 0, 

At;V„-— — m(a„)— =0 (27) 

/ ICJjn elVl 

We can re-write equation (!27|) in terms of the neutrino mass m{am) and dark energy density 
([5]) as, 

-log(^(l-a)=^i^ (28) 



Equation (128|) has the solution 



and 

Tm(cr„ 

The neutrino mass fl24l) reduces to the form, 



(29) 



(30) 



Mc'm) = j^e (31) 

Substituting for m((Tm) in (|30|1 we can solve for the mass parameter M, 

M = = 3.0 X lO^^GeV (32) 

V '2 Pa 

Substituting the expression ( 132|1 for M in equation ( l29l) we can solve for e which turns out 
to be 

1.0 X 10"^ (33) 



Using 032]) and f l5^ in the the expression fl?T]) . the neutrino mass turns out to be, 

eV 

m(a^) = 2.0 X 10"^^ (34) 
V A 

which is consistent with the assumption that atleast one species neutrinos is relativistic in 
the present era. 



4.2. Relativistic non-decoupled neutrinos 

Now we consider the case of neutrinos which are still in thermal equilibrium with some 
massless hidden sector radiation and which are still relativistic. The thermodynamic po- 
tential of relativistic neutrinos in terms of neutrino mass m(0) is obtain from ( llOp in the 
relativistic limit m < T is 



(27r)3 3 V(p2 + m2)V2y exp((p2 + m2)i/2/r) + 1 

and expanding the integrand in powers of m/p. Taking Qi = 2 for one relativistic neutrino 
species we obtain in the leading order in m/T, 

^.(</') = -^. = -^2^' + ^^^ (36) 
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Then the total thermodynamic potential is 

= Xv^a^ + —miaf - —T^ (37) 
^ ' 24 ^ ^ 360 ^ ^ 

The minima of Q occurs at the value of a = am given by the equation Q' = 0, 

A.V„-^mK)^=0 (38) 

We can re-write equation ( 138|1 in terms of the neutrino mass m{am) and dark energy density 
© as, 



.og(r^(l-.^)J^^= (39) 



Equation (l39l) has the solution 



and 

Tm(cr„ 



The neutrino mass (I23D is of the form 



(40) 



(41) 



2 

M 



Substituting (H2!) for m{am) in fj4T|) we can solve for the mass parameter M, 

M = , = 5.2 X lO^^GeV (43) 

Substituting the expression (H3|) for M in equation (HOj) we can solve for e which turns out 
to be 

1 

3.0 X 10"^ (44) 



24 

Using ( H3|) and dH]) in the expression P2|) we find that the neutrino mass is 

eV 

m{am) = 6.9 x 10"^^ (45) 
V A 

This is smaller than the neutrino temperature T = 1.710~^eV which is consistent with the 
assumption that the neutrino species is relativistic at the present epoch. 
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5. CONCLUSIONS 



In this paper we have studied the possibihty that observed the dark energy is due the 
shift of the Higgs field from the zero-minima of the Higgs potential due to the pressure 
of neutrinos. We have studied four different scenarios for the cosmological neutrino back- 
ground. Neutrinos could be relativistic and decoupled, non-relativistic and decoupled as in 
the standard cosmology. We have also considered the possibility that neutrinos are still in 
thermal equilibrium with some hidden sector radiation 2l(| for comparison. For our model 
to work the neutrino mass as function of the Higgs field must have a very steep gradient. We 
have considered logarithmic neutrino mass functions which meet this criterion. Interestingly 
the same mass model works for relativistic neutrinos whether they are decoupled or inter- 
acting. The same is true for the non-relativistic case, one neutrino mass model works for 
both the decoupled and the interacting neutrino scenarios even though the non-relativistic 
interacting neutrinos have a Boltzmann suppression in the number density. The neutrino 
mass function has a mass scale M, and dimensionless small parameter e as free parameters. 
For the case of non-relativistic neutrinos M = 3 x lO^^GeVA^/^ while for the case of relativis- 
tic neutrinos M ~ lO^^GeV. The GUT scale comes out naturally owing to the numerical 
relation {p\Y^^ ^ jM. Neutrino which is consistent with this scenario turns out to be 
nil, = 2 X 10"'^ eV (for the non-relativistic cases). This is in the mass range consistent with 
the solution of the atmospheric neutrino problem.The parameter e ~ 10~^ — 10~^ is put in 
to ensure that the minima-of the Higgs-neutrino thermodynamic potential is achieved at the 
present era. 

In this paper we explain the dark energy in terms of known particles. The neutrino mass 
must be a steep function of the Higgs field in order achieve this. It would be interesting to 
see if there are well motivated particle physics theory which gives rise to such neutrino mass 
functions. 
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